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Abstract 

Crystals are paradigms of ordered structures. While order was once seen as synony¬ 
mous with lattice periodic arrangements, the discoveries of incommensurate crystals and 
quasicrystals led to a more general perception of crystalline order, encompassing both 
periodic and aperiodic crystals. The current definition of crystals rests on their essen¬ 
tially point-like diffraction. Considering a number of recently investigated toy systems, 
with particular emphasis on non-crystalline ordered structures, the limits of the current 
definition are explored. 


1 What is order? 

The human brain is very skilled at detecting patterns and recognising order in a structure, and 
ordered structures permeate cultural achievements of human civilisations, be it in the arts, 
architecture or music. The ability to detect and describe patterns is also at the basis of all 
scientific enquiry; see Mumford & Desolneux (2010) for more on pattern theory. It may thus 
be surprising that a concept as fundamental as order does not have any well-defined precise 
meaning, and that it appears to be rather challenging to come up with a proper definition 
of what constitutes order in a structure. As a consequence, there currently is no satisfactory 
measure to quantify order in any given spatial structure. 

There are two common approaches to tackle this issue. One is to employ diffraction, which 
effectively measures two-point correlations in the structure; see Cowley (1995) for background. 
For kinematic diffraction, in the far-held approximation, the diffraction measure is the Fourier 
transform of the autocorrelation (or Patterson) function. Diffraction is the approach taken 
to characterise crystalline materials. The current dehnition of a crystal, which is based on 
its diffraction, derives from a dehnition which hrst appeared in the terms of reference of the 
lUCr Commission on Aperiodic Crystals, published in the 1991 report of the lUCr Executive 
Committee (International Union of Crystallography, 1992, p. 928); in fact, it is less general 
than what the commission proposed. The following quotes the more specihc dehnition given 
in Authier and Chapuis (2014), and used in the lUCr Online Dictionary of Crystallography. 

A material is a crystal if it has essentially a sharp diffraction pattern. The word 

essentially means that most of the intensity of the diffraction is concentrated in 
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relatively sharp Bragg peaks, besides the always present diffuse scattering. In 
all cases, the positions of the diffraction peaks can be expressed by 

n 

in >3). (1) 

i=l 

Here a* and are the basis vectors of the reciprocal lattice and integer coefficients 
respectively and the number n is the minimum for which the positions of the peaks 
can be described with integer coefficient h^. 

The conventional crystals are a special class, though very large, for which n = 3. 

The prominent role of the word essentially shows that this is a humble definition, in the 
sense that it reflects our limited knowledge of the structures we may potentially encounter 
in nature. The interpretation given in the definition that ‘essentially’ means that most of 
the intensity is concentrated in Bragg peaks means that the integrated contribution from the 
background must be weak compared to the Bragg diffraction, which is rather arbitrary, as 
any Bragg diffraction indicates order. By allowing the integer n to be larger than the three 
space dimensions we live in, aperiodic crystals are included in this definition, and conventional 
(periodic) crystals have become a special class (for which n = 3). Note that n is restricted to 
be finite here, so this particular form of the definition excludes pure point diffractive systems 
with non-finitely generated Fourier modules (over integer coefficients); the definition stipulates 
that Bragg peaks in crystals can be indexed by a finite number of integer coefficients. Note 
that the definition originally proposed in 1991 did not include this restriction (International 
Union of Crystallography, 1992, p. 928). 

Because the inverse problem of diffraction is inherently difficult (Bombieri &: Taylor, 1986) 
and, in general, not unique (Patterson, 1944), we do not have a complete characterisation 
of structures that show pure point diffraction (which means that the diffraction consists of 
Bragg peaks only), even in the idealised case of an ideal, perfect structure. Neither do we 
have a good understanding of what structures with an essentially pure point spectrum may 
look like. 

The second approach, which is particularly suited to stochastic systems, employs the 
entropy of a structure. Entropy takes into account the number of different local configurations 
of a system, and how this number grows with the system size; normally you are looking at an 
exponential growth with the system size, and any sub-exponential growth corresponds to zero 
entropy. Clearly, entropy can distinguish deterministic from random systems, and looking at 
different forms of scaling behaviour makes it possible to differentiate, at least to some extent, 
between different degrees of disorder. However, any deterministic structure has zero entropy 
(as has any sufficiently small deviation from it), so entropy is a rather crude measure of order. 

In this article, the current state of knowledge of mathematical diffraction of structures 
is summarised and discussed in relation with our notion of crystalline order. The current 
article attempts to present a broad overview only; for details on calculations and further 
background, we refer to recent survey articles (Baake & Grimm 2011a, Baake & Grimm 
2012) and to the monograph by Baake & Grimm (2013), as well as to references therein. 
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Using a number of explicit example structures with different types of diffraction spectrum, 
the range of possibilities is explored, contributing to the ongoing discussion on what order 
means in crystals and beyond (van Enter & Mi§kisz 1992, Lifshitz 2003, 2007, 2011). 

2 Mathematical diffraction 

In 2014, we were celebrating the international year of crystallography, and were looking back 
at a century of rapid developments in crystallography since von Laue (Friedrich, Knipping &: 
von Laue 1912, von Laue 1912) and father and son Bragg (Bragg & Bragg 1913) first employed 
X-ray diffraction to analyse the atomic structure of crystalline materials; see Authier (2013) 
for a historical account. In the simplest setting, which is suitable in particular for X-ray 
diffraction, it is sufficient to describe the kinematic scattering of radiation by the sample, 
and consider the far-held (Fraunhofer) limit for the outgoing radiation. The calculation of 
the diffraction pattern of a given structure then becomes possible by means of harmonic 
analysis, while the corresponding inverse problem of determining a structure from its pattern 
of diffraction intensities is, in general, difficult and non-unique, even in this simplihed setting. 
This section attempts to present a summary of mathematical diffraction theory, highlighting 
the ideas and the havour of the approach without going into technicalities, while trying to 
explain some of the technical terms by means of simple examples and familiar notions; for 
mathematical details, the reader is referred to Baake & Grimm (2013). 

2.1 What is a measure? 

A mathematically satisfactory approach to describe extended (infinite) systems, such as ideal 
crystals, is provided by using measures to describe both the distribution of matter in the 
scattering medium and the distribution of scattered radiation intensity in space. In math¬ 
ematics, measures are the natural concept to quantify spatial distributions, and are related 
to the notion of integration. The general approach to measures in mathematics is rather 
technical, but there is a simpler way to think of measures (which is due to a result called 
the Riesz-Markov representation theorem; see Reed & Simon (1980) for details). Indeed, it 
is possible to regard a measure as a linear functional, which is a linear map that associates a 
number to each function from an appropriate space. A familiar example is the integral of a 
function, which is the example we start with. 

A well-known and widely used measure in mathematics is Lebesgue measure, which is 
commonly used in integration of functions on the real numbers M. We denote Lebesgue 
measure by the letter A. If /: x f{x) is a function on R, the Lebesgue measure of / is 

Hf) = [ f{x)dX{x) = [ fix)dx, 

where the usual shorthand dx is used for integration with respect to Lebesgue measure. So 
Lebesgue measure associates to a function / a number, which is its integral. 

The Lebesgue measure of a set A C R, written as A (A), is given by 

A(A) = A(ly 4 ) = f lA{x)dx = f dx = vol(A), 

Jm. J A 
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where is the characteristic function (or indicator function) of A, which takes the value 
l^(x) = 1 for all X G A, and l^(x) = 0 otherwise. The Lebesgue measure of a set is what we 
call the volume (as in this case we are in one dimension, the volume is a length); for instance, 
the Lebesgue measure of the interval / = [a, b] with 6 > a is X{I) = A([a, b]) = b — a. Lebesgue 
measure is the unique translation-invariant measure on M (meaning that A(yl) = A(^ -|- 1) for 
any translation t E M) that assigns the volume 1 to the unit interval [0,1]. Lebesgue measure 
on M generalises in the familiar way to d-dimensional space corresponding to d-dimensional 
(multiple) integrals. For simplicity, we will mainly consider the case d = 1 in what follows. 

Another well-known and important measure is the Dirac measure (or point measure) at a 
point X G M, denoted by 5x- It describes a localised structure at a point x in space, with total 
measure 1. This means that, if / is a function on R, its point measure at x is Sx{f) = fix). 
In the physics literature, the point measure is often written like a function 5{x) (which can be 
considered a generalised function or distribution obtained as a limit of functions, for instance 
of a sequence of Gaussian functions of integral 1, centred at the origin, and with a decreasing 
width, which then become increasingly sharper), with the suggestive notation 

4(/) = / fiy)Six-y)dy = /(x). 

Jr 

This notation can be used consistently as long as one remembers that Dirac’s 5 is not a 
function in the usual sense. As above, one can also define the point measure of a set A C R, 
which is 5xiA) = 5(1^) = ^Aix), so 5xiA) = 1 if x G A and 0 otherwise. 

2.2 Dirac combs 

Point measures are often used to describe a set of localised scatterers in space. Given a set of 
scatterers located at points in a subset A C R (which we usually assume to be a Delone set, 
which means that it neither contains points that are arbitrarily close to each other nor holes 
of arbitrary size), we can associate a measure 

< 5 ^ := 4 

x&A 

which we call the Dirac comb (a term coined by de Bruijn (1986); see also Gordoba (1989)) 
of A. An example of such a Dirac comb is which is the uniform Dirac comb 

on the integer lattice. 

By introducing scattering weights rc(x) at position x G R (which in general can be complex 
numbers, but we will assume to be real for the purpose of this exposition), we arrive at a 
weighted Dirac comb 

uj = w 5 a = Y^w(x)5x, ( 2 ) 

xeA 

which can serve as a model representing a scattering medium containing different types of 
scatterers. Any measure of this type, consisting of a (weighted) sum of point measures, is 
called a pure point measure (with respect to Lebesgue measure). It is possible to include 
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Figure 1: Schematic representation of the weighted (periodic) Dirac comb a; of Eq. ([3|). Point 
measures a 5x are represented as columns at position x of height proportional to their weight 
a. 

realistic scattering profiles by considering convolutions with appropriate motives, so this ap¬ 
proach is not as restrictive as it may seem. A schematic representation of an example, the 
weighted (periodic) Dirac comb 

^ = 57 . +\ + i (3) 

is shown in Figured) 

Attaching scattering profiles to a Dirac comb is one way to represent a continuous scat¬ 
tering intensity in space. Of course, there is a more direct approach if the scattering intensity 
is described by a continuous distribution is space. If ^ is a such a continuous (or, at least, 
locally integrable) function on M, it dehnes a measure /U on M via 

= [ f{x)dn{x) = [ f{x)Q{x)dx. 

In this case g is called the density of the measure //. Any measure g that can be written in 
this form is called an absolutely continuous measure (with respect to Lebesgue measure). 

The measures we are interested in are those which describe distributions (of scatterers or 
scattering intensity) in space, and one physical restriction we would like to impose is that any 
finite region of space can only contain a finite total scattering strength or finite intensity. The 
measures that satisfy this property are called translation bounded measures. A Dirac comb 
8a of a Delone set A C M is always translation bounded, because there can be only finitely 
many points in any finite region of space, due to the minimum distance between points. The 
same is true for a weighted Dirac comb provided that the weight function w{x) is bounded. 
An examples of a measure that is not translation bounded would be a Dirac comb of a point 
set with an accumulation point, for instance X]nez\{o} ^i/n- measure, any interval 

containing the origin contains infinitely many point measures, and thus has infinite measure. 

2.3 Lebesgue decomposition 

A central result in measure theory is Lebesgue’s decomposition theorem. It states that any 
measure can be written as a sum of three components in a unique way (with respect to a 
reference measure, which is our case will always be Lebesgue measure). If ^ is a measure on 
M, the three components are denoted as 


A — h'PP /^sc T hac 
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Figure 2: Sketch of the classic middle thirds Cantor set construction (inset) and the ‘Devil’s 
staircase’ distribution function F of the corresponding singular continuous measure. 


and are called the pure point component /Xpp, the singular continuous component psc and 
the absolutely continuous component pac- We have met typical examples of pure point and 
absolutely continuous measures above, so the obvious question is what a singular continuous 
measure looks like. As it is defined, it is all that is ‘left’ if you remove the pure point part 
(consisting of a sum of weighted point measnres) and the absolutely continuous part (which is 
described by a locally integrable density function) — but this does not really help to gain an 
understanding of what such a measure represents. Singular continuous measures are rather 
weird objects indeed; they do not give weight to any single point in space (because otherwise 
it would have a pure point component), but are concentrated on sets of vanishing volume 
(because otherwise you could describe part of it by a density function). 

The best-known example of a singular continuous measure is provided by the classic 
middle-thirds Cantor set; see Figure [2j Starting from the unit interval, the Cantor set is 
constructed by removing the middle third of it, then removing the middle thirds of the two 
resulting intervals, and iterating this procedure ad infinitum. The corresponding Cantor 
measure is constructed by starting from the Lebesgue measure on the interval, so we have 
total measure 1, and at each step distributing the mass equally onto the constituent intervals. 
In the limit, the total measure is thus still 1, but there is neither any isolated point that 
carries a finite measure (since the measure of each interval in the nth step is 2“”, so it 
vanishes in the limit) nor any interval of finite length that is in the support of the measure 
(meaning that the measure does not vanish on it). The measure constructed in this way is 
thus singularly continuous, and can be described in terms of a distribution function which is 
a ‘Devil’s staircase’. This function is constant almost everywhere, and displays a hierarchy 
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of plateaux in its graph (see Figure [2]) which reflect the hierarchy of gaps produced by the 
excision steps in the Cantor construction. This function plays the role of the integrated 
density for the singular continuous measure. 

Lebesgue’s decomposition provides a rigorous way to separate the diffraction measure 
of a structure into its pure point (Bragg) part and its singular and absolutely continuous 
components. However, using this as the definition really only makes sense if one works with 
infinite systems (because finite systems will always have absolutely continuous diffraction). 
This is similar in spirit to the definition of a phase transition in materials (as a discontinuity 
in a thermodynamic potential), which again only applies in the mathematically rigorous 
sense to an infinite system (because for finite systems these potentials are smooth functions). 
Nevertheless, these concepts have proved useful for applications to macroscopically large 
(albeit finite) systems. 


2.4 Autocorrelation and diffraction 

A key quantity in diffraction theory is the autocorrelation, which quantifies the two-point cor¬ 
relation of a structure. In crystallography, this is often called the Patterson function. If the 
material is described by a (real) density function g on R (so we deal with an absolutely con¬ 
tinuous distribution), the autocorrelation is an absolutely continuous measure whose density 
is the familiar convolution 

P{x) = [ g{y)g{y + x)dx = {g*g){x), 

Jm. 


where g is the function defined by g{x) = g{—x). 

In the case that the material is described by a one-dimensional weighted Dirac comb u 
of the form given in Eq. ([2]) (with real weight function w{x)), the autocorrelation is a pure 
point measure 

zeA-A 


with non-vanishing contributions only at distances z in the difference set A — A = {x — y \ 
X, 7 / G 4} (which you may interpret as the set of interatomic distances). The point masses for 
interatomic distances z are weighted by the autocorrelation coefficients yiz), which are given 
by 


r]{z) 


lim —— 
ij—>oo 2i? 


w{x) w{z 

\x\<-R.,x^A 

z—x^A 


provided that these limits exist. Note that 2R is the length of the interval {—R,R), so the 
autocorrelation coefficient r]{z) is precisely the volume-averaged two-point correlation for the 
interatomic distance 2 ;. 

Using the language of measures, these equations can be neatly summarised as follows. 
Given a (translation bounded) real measure cj, again for simplicity in one dimension, its 
autocorrelation measure 7 is defined as 


'J = LJ ® CO 


lim 

R^oo 


^\r * ^\r 

2R 


(4) 
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provided the limit exists. Here, denotes the restriction of a; to the interval {—R,R), and 
J1 is defined via Jl{g) = fi{g) with g{x) = g{—x) as above. Finally, * denotes the standard 
convolution of measures, which is the appropriate generalisation of the convolution of func¬ 
tions. The autocorrelation of u is thus the volume-averaged convolution @ (also called the 
Eberlein convolution) of uj with its ‘hipped-over’ version a;, and thus picks out the two-point 
correlations in the structure described by uj. This approach to mathematical diffraction was 
pioneered by Hof (1995). 

The diffraction measure is then the Fourier transform 7 of the autocorrelation, so essen¬ 
tially it provides a spectral analysis for the two-point correlations in the original structure. 
It is a translation bounded, positive measure, which quantihes how much of the kinematic 
scattering intensity reaches a given volume in space. Lebesgue’s decomposition 

T Tpp T Tsc ~k Tac 

into its pure point part (the Bragg peaks, of which there are at most countably many), its 
absolutely continuous part (the diffuse background scattering, described by a locally integrable 
density function) and its singular continuous part (which encompasses anything that remains) 
provides a mathematically rigorous definition of the different types of diffraction. For the 
definition of a crystal cited above, it is the pure point part 7 pp that matters — a crystal is a 
structure where this part represents the majority of the diffracted intensity (there will always 
be some continuous diffraction in practice), though this alone does not guarantee that the 
positions of Bragg peaks can be indexed by a finite number of integers. Indeed, there are 
examples of systems that are pure point diffractive, meaning that 7 = 7 pp, where this is not 
the case; we shall meet an example below. 

3 Periodic crystals 

A conventional, periodic crystal is described as a lattice-periodic structure, corresponding to 
an ideal infinite crystal without defects or surfaces. A periodic crystal is characterised by its 
periods (translations that keep the crystal invariant), which form a lattice R (because any 
linear combinations of periods are also periods), and by the distribution of scatterers in a unit 
cell (fundamental domain) of this lattice. Here, a lattice F in d-dimensional spac^ is the set 
of integer linear combinations of d linearly independent basis vectors, so it can be written in 
the form 

^ ■= I Uj G z| , 

where vi G for 1 < f < d, are linearly independent vectors in W^. Familiar examples are 
the integer lattice Z in one dimension, the square lattice in two dimensions and the simple 
(primitive) cubic lattice Z^ in three dimensions. 

If the scattering medium has a (periodic) crystal structure described by a lattice T, it can 
always be represented as a measure 

u = fj, * 6p , (5) 

^Because the generalisation to higher dimensions is straightforward, we give the results for the general case, 
although you can always think of cases with d < 3; see also the examples below. 




where ^ can be chosen as a finite measure which describes the decoration of the fundamental 
domain, while the Dirac comb 5p ensures lattice periodicity. 

The corresponding autocorrelation measure is a T-periodic measure that can be calculated 
from the appropriate generalisation of Eq. Q as 


7 = dens(T') (n* fi) * 6 p , 


( 6 ) 


using 6 p = 6 p and 6 p® 6 p = dens(T) 6 p, where dens(T') denotes the density (per unit volume) 
of the lattice F. The diffraction measure 7 is then given 


7 = (dens(T))' 


(7) 


This provides the familiar result for periodic crystals: Any perfect periodic crystal with lattice 
of periods F shows pure point diffraction with Bragg peaks located on the reciprocal lattic^ 
F*, and the intensity of the Bragg peak is determined by the density of the crystal lattice 
F and by the continuous function ml , which depends on the distribution of scatterers in a 
fundamental domain of F. By expressing the reciprocal lattice positions as linear combinations 
of basis vectors of the reciprocal lattice F*, this can be cast in the form of Eq. © with n = d. 

As a one-dimensional example, consider the weighted Dirac comb w of Eq. ([3|) and Eigure[T] 
It can be written as 


UJ = 


+ 4 + 2 4 


2 <53) * Sz , 


so here a = 60 + \5i + i(5i -|- 4(^3 describes the four scatterers within the fundamental domain 
^4 ^2 ^4 

[0,1) of the lattice T = Z. Note that in this case we have 


OJ = 


T 4 i T 2 1 + 4 3 

^ 4 ^ 2 ^ 4 


= OJ . 


due to the equivalence of positions differing by integers in the periodic lattice and the sym¬ 
metric distributions of scatterers in the fundamental domain. Let us now calculate the auto¬ 
correlation and diffraction of this comb. 

Clearly, since all distances are multiples of the autocorrelation in this case will have a 
similar form as the comb w itself, just with different coefficients, which are given by summing 
up the products of the weights of scatterers with a given separation. To obtain these coeffi¬ 
cients, one can compute the convolution of oj with itself (or equivalently of ^ with itself) term 
by term, using the relation 5x* dy = dx+y This gives 


7 = (t* + Ti + !7 + i7)*'5z 


Eor instance, the coefficient ^ -|- (i)^ -|- (i)^ -|- (i)^ of Jq comes from adding up the 

contributions to integer distances. A schematic presentation of the autocorrelation 7 is shown 
in Eigure [3l 


^This follows from Eq. © by an application of Poisson’s famous summation formula, which can be cast 
as 5p = dens(T) Sp *, where F* denotes the dual (or reciprocal) lattice of F ; see Baake & Grimm (2013) for 
details. 

®Note that we prefer to define the Fourier transform of a function / as f{k) = f{x) Ax. Due to the 

factor 277 in the exponent, one avoids the appearance of such factors in the definition of the reciprocal lattice. 
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Figure 3 : Schematic representation of the autocorrelation (top) and diffraction (bottom) of 
the weighted Dirac comb w of Eq. ([ 3 ]) and Figure [H 

The corresponding diffraction measure 7 is obtained by taking the Fourier transform, 
using that a6x = This gives 

7 = (^ + |(-l)" + icos(f))*5z 

= 4 64Z + < 54 Z +2 + 3 <^4Z+{1,3} • 

The diffraction pattern is thus periodic, but with period 4 (due to the smallest distance 
between scatterers being 4 ). A schematic picture of the diffraction pattern is shown in 
Figure El 

As a second example, consider a two-dimensional crystal with lattice of periods F = 1 ?, 
with two scatterers (of unit scattering strength) per unit cell, one placed at lattice points 
and the other at an arbitrary position (a, 6) € [ 0 , 1 )^- The corresponding point set is A = 
U (7? + (a, 6)), and the Dirac comb can be written as a; = * (5^2 = (Jq.o + ^a,b) * <5z2- Its 

autocorrelation is 

7 = {q*q) * 

= (<^(0,0) + ^{a,b)) * (<^(0,0) + ^-{a,b)) * 

= (2<I(0,0) + ^{a,b) + ^-{a,b)) * <^22 • 

The corresponding diffraction measure is then 

7 = \Q\'^{kJ)Sz2 

= (2 + 2Re(e-2"(^“+^^)))5^2 
= (2 -|-2 cos(( 27 r(/ca + £6))) 6^2 

= (2 cos [it {ka + 6.^2 

for k,£ £ Z. Note that, while the diffraction measure is supported on as expected (as 
Z^ is self-dual), it need not have any non-trivial period. In fact, 7 is periodic in one or two 
directions precisely if one or both coordinates a and b are rational, respectively; an example 
with one periodic direction is shown in Figure IH Although the positions of Bragg spots for a 
lattice-periodic structure are again lattice-periodic, in general the intensities will not respect 
the periodicity of the dual lattice. 
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Figure 4: The left panel shows a schematic representation of the two-dimensional toy crystal 
discnssed in the text, with two scatterers of equal strength at positions ( 0 , 0 ) and (a, b) = 
(ij ■^) of the fundamental domain (indicated by shading). The corresponding autocorrelation 
7 is shown in the central panel, while the corresponding diffraction measure 7 is displayed 
in the right panel. Here, a point measure is represented by a dot that is centred at the 
position of the peak and has an area proportional to the weight of the point measure. The 
irrational position in the vertical direction leads to an incommensurate modulation of the peak 
intensities in this direction, which the diffraction is periodic with period 3 in the horizontal 
direction. 

4 Aperiodic crystals 

The arguably best understood class of aperiodic structures are cut and project sets, also called 
model sets. Model sets can be viewed as a natural generalisation of the notion of quasiperiodic 
functions, which goes back to Harald Bohr (1947), and were first introduced by Yves Meyer 
(1972) in the context of harmonic analysis. The basic idea of the construction is to obtain an 
aperiodic structure as a suitable ‘slice’ of a higher-dimensional periodic lattice, which is then 
projected onto a suitable space of the desired dimension. For simplicity, we mainly consider 
the case where the higher-dimensional space is a Euclidean space of the form with 

being the physical space (sometimes also called the direct or the parallel space) that hosts 
the aperiodic structure (so 1 < d < 3 for physically relevant cases), and the internal (or 
perpendicular) space which is used in the construction. 

Let us start with an example, where d = m = 1. In this case, we are projecting a one¬ 
dimensional aperiodic structure from a two-dimensional (periodic) lattice. In this example, 
the lattice L is given by all integer linear combinations of two basis vectors, which we choose 
as ( 1 , 1 ) and (r, 1 — r), where t = {1 + is the golden ratio, so we have 

£ = Z (1,1) + Z (t, 1 — t) = {m (1,1) -|- n (r, 1 — r) | m, n € Z} . 

The lattice points are shown as black dots in Figure [5j The lattice is oriented such that 
the horizontal space along the ( 1 , 0 ) direction is the physical space and the vertical direction 
along (0,1) corresponds to the internal space. We call the projection to the physical space vr, 
and the projection to the internal space vTint- The projections of all lattice points L = 7r(£) 
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to physical space and L* = 7rint(^) to internal spacqj are both dense in their corresponding 
one-dimensional spaces. The set L is explicitly given by L = Z[r] = {m + nr \ m,n ^ Z}, so 
all integer combinations of multiples of 1 and r, which is dense because r is irrational, and L* 
has the same form. Note that the projections are one-to-one in both directions. In particular, 
any point in L corresponds to a uniquely defined point in £. In fact, 7r“^(x) = {x,x*), where 
the *-map is defined by mapping 1 to 1 and r to I — r (which corresponds to the ‘algebraic 
conjugation’ \/5 —V^), so (m -|- nr)* = m -|- n(l — r) = m -|- n — nr for all m,n £ 1^. 

The final ingredient that we require is a ‘window’ W in the internal space, which we 
choose to be the half-open interval W = (—l,r — 1]. Shifting it along the physical space 
sweeps out the shaded horizontal strip in Figure [SI The lattice points that fall within this 
strip produce the set {x £ C \ 7rjj^^(x) £ W}, and their projection onto the physical space 
thus A = {7r(x) \ X £ C and 7rjj^^(x) £ W}. Using tt{C) = L and the *-map, this point set can 
equivalently be written as 

A = {x £ L\x* £W}. (8) 

Sets of this form are called cut and project sets or model sets. The condition that x* £W 
selects a discrete subset of the dense set point set L, in fact, a very special discrete subset 
where points are separated either by intervals of length 1 (for short intervals s) or by intervals 
of length r (for long intervals tj. As it turns out, this projection yields the famous Fibonacci 
sequence ... tsilstsi ... of long {€) and short (s) intervals, which can be generated by the 
two-letter substitution rule I ^ Is^ s ^ t. In particular, dividing the window into two parts 
as follows 

Ws = (r — 2,r—1] and Wi = (—1,t — 2] 

shows that the sets of left endpoints of short or long intervals are given by the projection of 
lattice points that fall within the corresponding sub-strip, so 

As = {x £ L \ x'^ £ Ws} and Ag = {x £ L \ x* £ , 

with A = AsA Af^. Hence the set of left endpoints of short or of long intervals separately are 
model sets with windows Wg and Wi, respectively, while the set of all left interval endpoints 
is a model set with window W = Ws U Wi\ compare Figure [5l 

As mentioned above, this construction can be generalised to physical and internal spaces 
of any dimension. The general cut and project scheme (CPS) for Euclidean model sets can 
be summarised in the following diagramme 

U U U dense 

7r(£) ^ (9) 

L -^^ L* 

Here, £ C is a lattice in the (d-l-m)-dimensional space x R”* = and vr and 

denote the natural projections from this space onto the physical and internal spaces R'’* and 

^Note the difference between the star symbol * used here and the * used for the dual or reciprocal lattice. 
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Figure 5: Schematic representation of a natural projection approach for the Fibonacci chain 
from the planar lattice spanned by the vectors (1,1) and (r, 1 — r). 

M”*, respectively. We assume that the point set L = 'k{C) C which is the projection of 
the lattice points into the physical space, is a bijective image of which means that no two 
lattice points in £ project onto the same point in L. In other words, each point in L can be 
‘lifted’ to a unique lattice point in £, and the inverse map 7r“^ is well-defined on all elements 
of L. This ensures that the star-map *: x x* is well-defined on L, so each point in L has 
a unique associate in internal space; see Moody (2000) for details. Finally, we assume that 
the corresponding set L* = 7r.^^(£) C in internal space is dense. 

Given a CPS, the second ingredient in the definition of a cut and project set is the window 
(sometimes also called an acceptance domain) W C M™, which is assumed to be a sufficiently 
nice subset of (technically, a relatively compact subset with non-empty interior). A cut 
and project set is then defined by selecting all points x in the projected lattice L whose 
companion x* in internal space falls inside the window W. Expressed as an equation, this 
means that any set of the form 


A = {x E L I X* E W}, (10) 

or indeed any translate of such a set, is what we call a model set. The technical conditions on 
the window W ensure that T is a Meyer set (Meyer 1972, Moody 2000), which means that 
the difference set 

A — A := {x — y \ x,y & A} 

is uniformly discrete (so different distances between point in the structure differ by at least a 
fixed amount) and that the set A is relatively dense (which means that there are no arbitrarily 
large ‘holes’ in the point set). If the boundary dW of the window W is nice in the sense that 
it has zero volume (in the sense of Lebesgue measure), we refer to A as a regular model set. 
The setting of Eq. ([9]) can be generalised further to allow for the internal space to be a locally 
compact Abelian group (Meyer 1972, Moody 2000, Schlottmann 2000). 
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Figure 6: Equivalent description of the Fibonacci chain in terms of ‘targets’, often referred 
to as ‘occupation domains’ or ‘atomic surfaces’. 

It is worth mentioning that there are various equivalent ways of interpreting the cut and 
project construction. One commonly used approach attaches an inverted copy of the window 
as a ‘target’ to each lattice point, and projected points are then obtained as the intersection 
of these targets with the physical space; see Figure E] for an illustration of the Fibonacci case. 
Albeit equivalent, this description offers a simpler way of interpreting experimental data, 
and is therefore the preferred presentation of the cut and project approach in experimental 
research papers, where it is often referred to as the section method. Apart from providing 
an intuitive meaning for the targets as ‘atomic surfaces’, this approach allows for additional 
variation (by deformations of the targets) that can be exploited, for instance in the description 
of modulated phases. For further variants of the cut and project method, the reader is referred 
to Chapter 7.5 in Baake & Grimm (2013). 

Familiar examples of model sets are some one-dimensional substitution-based structures 
such as the Fibonacci chain discussed above, and some of its generalisations. Planar examples 
include the Penrose tiling and the Tubingen triangle tiling with decagonal symmetry, the 
Ammann-Beenker tiling with octagonal symmetry and the shield tiling with dodecagonal 
symmetry, which can all be obtained by projection from four-dimensional lattices. Structure 
models of icosahedral quasicrystals usually employ model sets based on the (primitive, face- 
centred or body-centred) hypercubic lattice in six dimensions. These examples have direct 
application to the crystallography of quasicrystals, and serve as models for the structure 
of decagonal, octagonal, dodecagonal and icosahedral quasicrystals, respectively; compare 
Steurer &: Deloudi (2009) for details. Realisations of model sets with other symmetries, such 
as planar sevenfold symmetry for instance, have as yet not been observed in nature, and the 
same is true for model sets where the internal space is not Euclidean. Nevertheless, such 
systems share many properties with the familiar quasicrystalline cases, and should not be 
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excluded a priori. Even if it may prove impossible to realise such structures in self-assembled 
systems, we may be able to produce these in purpose-made manufactured structures at various 
length scales, from macroscopic to nanometre and atomic scales. 

Arguably the most important result in the theory of model sets, in the context of math¬ 
ematical diffraction theory, is the proof that regular model sets have pure point diffraction. 
Three different approaches have been used to prove this statement. The first proof using 
methods of dynamical systems theory was completed by Schlottmann (2000), employing an 
argument by Dworkin (1993) and the mathematical diffraction approach of Hof (1995); see 
also Lenz & Strungaru (2009) for further developments. An alternative approach employs 
the theory of almost periodic measures (Baake & Moody 2004, Moody & Strungaru 2004, 
Strungaru 2005). Following a suggestion by Lagarias, Baake &: Grimm (2013) present a proof 
based on Poisson’s summation formula for the embedding lattice in conjunction with Weyl’s 
lemma on uniform distribution, which exploits the uniform distribution of projected lattice 
points in internal space. Although it has not been developed into a proof, there is also a 
somewhat complementary approach based on an average periodic structure; we refer to the 
recent review by Wolny, Kozakowski, Kuczera, Strzalka & Wnek, A. (2011) and references 
therein for details. 

Essentially, the pure point diffraction of a model set is a consequence of the underlying 
higher-dimensional lattice periodicity. Let us first discuss the example of the Fibonacci model 
set A of Eq. ([ 8 ]); compare also Figure [3 The pure point diffraction pattern is obtained again 
as a projection to physical space, but this time of the dual (or reciprocal) higher-dimensional 
lattice C*. In our case, this is the lattice generated by the dual basis vectors (r — 1 ,t) 
and ^^^(1, —1). The corresponding Fourier module is then 

L® = 7r(r) = ^Z[r], 

where Z[r] = {m -|- nr | m, n G Z} as above. The determines the positions of Bragg peaks, 
but what about their intensities? It turns out that the intensity is a function of the distance 
of the projected lattice point from the physical space, and roughly the larger the internal 
coordinate the smaller the intensity. The function in question is the absolute square of the 
Fourier transform of the window function (the characteristic function of the window), which 
is the function that takes the value 1 on the window and 0 otherwise. Its Fourier transform 
is 

A(k) = ^ sinc(7rTfc*), ( 11 ) 

5 

where sinc(x) = sin(x)/x, and k* is the image of k under the ★-map introduced above. A 
sketch of the diffraction pattern is shown in Figure [71 

Let us now return to the general result. For a regular model set A with Dirac comb 5^, 
the diffraction measure 7 can be written as 

7 = J;|A(A:)|2 4 , (12) 

k&L® 

where L® = 7 r(i 2 *), the projection of the higher-dimensional dual lattice, is the corresponding 
Fourier module on which the pure point diffraction is supported. For a Euclidean model set 
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Figure 7: Sketch of the projection of the dual lattice points giving rise to Bragg peaks in 
the diffraction pattern for the Fibonacci point set A, with scatterers of unit weight at all 
points. The function displayed on the right-hand side is the intensity function |^(A:)p. The 
Bragg peak at 0 has height (dens(yl))^ = (r + l)/5 0.5206, and the entire pattern (once all 

reflections are included) is reflection symmetric. 


with the CPS ([9]), £ is a lattice in a Euclidean space and the Fourier module L® is thus 

finitely generated, with rank d -h m. By choosing appropriate generating vectors, the pure 
point diffraction of Eq. (I12p can thus be recast in the form of Eq. ([1]) with n = d+m. However, 
in the general situation, where the internal space can be any locally compact Abelian group, 
this is not necessarily the case, as the Eourier module L® = vr(£*) does not have to be finitely 
generated. Note that the latter case is not covered by the definition of a crystal cited above, 
while it does include any model set based on a Euclidean CPS. 

The diffraction amplitudes A{k) are obtained by the Eourier transform of the characteristic 
function of the window W, 


m 


dens(yl) 

vol(kF) 






(13) 


According to Eq. (fT^ . it is the absolute square of these amplitudes that determine the 
intensity of a Bragg peak as position k G L®, with k* denoting the corresponding point in 
internal space. Eq. (I13|) gives the result for Euclidean model sets, the only difference for the 
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general case is that the volume (with respect to Lebesgue measure in Euclidean space) is 
replaced by the suitable invariant measure (Haar measure) on the locally compact Abelian 
group. 

5 Order beyond crystals 

The current definition of crystals thus covers conventional periodic crystals, incommensurate 
crystals and quasicrystals, and hence all currently known realisations of perfectly ordered 
structures in nature. While the question asked by Bombieri & Taylor (1986) has not yet been 
satisfactorily answered, it is clear that pure point diffraction is a severe constraint on the 
possible structure (Baake, Lenz &: Richard 1997), and recent results by Lenz &: Moody (2009, 
2011) indicate that model sets play a major role in a potential abstract parametrisation of the 
inverse problem. However, there are clearly well-ordered structures that do not possess this 
property, and this section will discuss a few characteristic examples. But first, we start with 
an example of a pure point diffractive system with non-finitely generated Fourier module, 
which thus possesses a diffraction pattern where Bragg peaks cannot be indexed by a finite 
number of integers. 

5.1 Pure point diffraction with non-finitely generated support 

Well-known examples of systems with non-Euclidean internal spaces are limit-periodic struc¬ 
tures. Let us explain this with the arguably simplest example, based on the period doubling 
substitution rule q \ 10,0 11, on the two-letter alphabet {0,1}. Any bi-infinite worcjfl 

w E {0,1}^ that satishes the hxed point property Q^iw) = w is specified completely by 
w{2n) = 1, w{4:n -|- 1) = 0 and w{4:n -|- 3) = w{n) for n E Z, while either letter can be chosen 
at position n = — 1. The two possible choices lead to two locally indistinguishable sequences 
(which means that any finite subsequence of one occurs in the other), and hence define the 
same system. 

The word w possesses a Toeplitz structure consisting of a hierarchy of scaled and shifted 
copies of Z which carry the same letter. Dehning the point set 

A = {re E Z I w{n) = 1} C Z (14) 

of the positions of the letter 1 in w, it is clear from the relations above that 2Z C A, as 
all letters at even positions are 1. But then, due to 'u;(4re -|- 3) = w{n), so are all letters 
w{8n+3) = w{2n) = 1, so 8Z-I-3 C A, and inductively one recognises that 2-4^Z-|-(4^ —1) C A 
for all integer i > 0. In fact, this hierarchy of scaled integer lattices describes the complete 
set, and we obtain the following representation for the set as a union (Baake, Moody & 
Schlottman 1998, Baake Sz Moody 2004, Baake & Grimm 2013) 

A = 2ZU(8Z + 3)U(32Z + 15)U... = |J ((2 • 4^ Z + (4^ - 1)) (15) 

e>o 

®Here and below, the notation A* denotes the set of bi-infinite sequences (..., a_ 2 , ci-i, Uq, bj, 02 ,...) with 
letters Ui, i chosen from a finite alphabet A. 
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Figure 8: Schematic representation of the Dirac comb 8a of the period doubling point set 
of Eq. (jl4h . All point measures have the same mass. The different shading highlights the 
Toeplitz structure, with point masses at even integers shown in black, point masses on 8Z + 3 
in dark grey and a single point mass in 32Z + 15 in lighter grey. 

of scaled (and shifted) lattices. Note that this result is for the case where we choose rc(—1) = 0 
(otherwise —1 has to be added to the right-hand side). A schematic representation of the 
corresponding Dirac comb 8a is shown in Figure [8l 

Using this representation for the point set A, the diffraction of the Dirac comb 8a can 
be computed explicitly; see Baake &: Grimm (2011b) for details. The scaled lattices with 
geometrically increasing period in the union in Eq. (|15l) give rise to Bragg peaks supported 
on the corresponding dual lattices, which are successively finer integer lattices scaled with 
the inverse factor. The diffraction spectrum is pure point, and the Eourier module can be 
parametrised as 


L® = Z[i] = I (r = 0,m G Z) or (r > l,m odd)} . 


(16) 


The diffraction measure is of the form of Eq. (I12jl , with diffraction amplitudes 


A /'m\ _ „2^ '’Trim 

or ® 


(17) 


at the positions in L®. The factor of | reflects the density of scatterers, as two thirds of 
positions are occupied. It is no coincidence that the model set expression applies — in fact, 
the set A can be described as a model set, but with a non-Euclidean internal space; in this 
case, the internal space is what is known as the space of 2-adic integers (which essentially 
consists of all fractions whose denominators are powers of 2, but with a specific definition of 


Figure 9: Schematic representation of the diffraction intensity pattern of the Dirac comb 
8a of Figure El The pattern is periodic with period 1 and consists of a dense set of Bragg 
peaks, where increasingly smaller peaks are located at rational numbers whose denominators 
are increasingly larger powers of 2. Note that only peaks corresponding to r = 0,1,2,3 in 
Eqs. (fTGll and (fTTll are visible here. 
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the distance of two numbers). A schematic representation of the diffraction pattern for the 
period doubling chain is shown in Figure [9l 

It is easy to generalise this example to other lattice-based substitutions in one or more 
dimensions; any substitution of constant length p with a coincidence in the sense of Dekking 
(1978), which means that the same letter appears at the same position for the images of 
all letters under a certain power of the substitution rule, is a good candidate, because it is 
always pure point diffractive and carries a natural p-adic structure. A well-known example of 
this type is the chair tiling, in its representation as a two-dimensional block substitution; see 
Robinson (1999) and Baake & Grimm (2013) for details. 


5.2 Order and singular continuous diffraction 

The paradigm of singular continuous diffraction is provided by the Thue-Morse system and 
its generalisations (Kakutani 1972, Baake &: Grimm 2008). Here, we consider a family of 
generalised Thue-Morse substitutions (Baake, Gahler &: Grimm 2012) 




J ^ p 1^ 


(18) 


on the two-letter alphabet {1,1}, where fe,!" G N and the case k = (. = 1 corresponds to the 
classic Thue-Morse case. Note that V denotes a string of j consecutive letters 1 here. The 
one-sided fixed point v = (y) satisfies 


'k^(k+tjm+r 


Vm, if 0 < r < k — 1, 

Vm, if k<r<k + i — 1 


(19) 


where m > 0 and 0<r<k + i — 1 and 1 = 1. It extends (by setting n_j_i = Vi for i > 0) 
to a symmetric bi-infinite fixed point word under the square of the substitution For 

instance, the symmetric bi-infinite fixed point for the classic Thue-Morse case k = i = 1 has 
core 

... Iiiliniliiniinnilinininiiliniliil... 

where the vertical bar denotes the origin. A schematic representation of the corresponding 
Dirac comb is shown in Figure [TOl 


> 


Figure 10: Schematic representation of the Dirac comb of the Thue-Morse chain with weights 
±1. Note that this is ‘balanced’ in the sense that positive and negative weights are equally 
frequent, so the average scattering strength is zero. 
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Figure 11: Approximating density functions (top) and corresponding approximating 

distribution function (bottom) for the diffraction of the classic Thue-Morse chain (k = 
i = 1) with N = 4 (left), N = 5 (centre) and N = 6 (right). 


The corresponding weighted Dirac comb on Z, interpreting the two letters as weights (with 
1 interpreted as —1), is thus given by a; = Its autocorrelation 7 = 

is also a Dirac comb on Z, where the autocorrelation coefficients rj{m) satisfy 7 ( 0 ) = 1, 
r]{—m) = r]{m) and the recursion relation 


r]{^{k + t)m + r) 


k + i 


(^a{k, £, r) r/(m)+ 
a{k, i,k + £ — r) rj{m + 1 )^ 


( 20 ) 


for arbitrary m € Z and 0 < r < k + i — 1. The recursion follows directly from Eq. (|19|) . 
with a{k, £,r) = k + £ — r — 2 min(/c, £,r,k -\- £ — r). This system has a unique solution, and 
properties of the solution show that the corresponding diffraction measure 7 is purely singular 
continuous; see Baake, Gahler & Grimm (2012) for the mathematical details of the argument. 

The diffraction measure is periodic with period 1 (due to the fact that the Dirac comb is 
supported on the integer lattice Z) and the diffraction intensity can be represented as a limit 
of a product, 

N 

= im{k+£rx), 

n=0 

which is known as a Riesz product, with 


??(x) 


1 + 


2 

k + £ 


k+l-l 

a{k,£,r) cos(27rrx). 

r=l 
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Figure 12: Continuous distribution functions for the diffraction of generalised Thue-Morse 
chains with i = 1 and k = 2 (left), k = 3 (centre) and k = 4 (right). 

The limit as N ^ oo has to be considered carefully. While the truncated product is a 
smooth function that can be interpreted as a density of an absolutely continuous measure, 
this is not the case in the limit, because it represents a purely singular continuous measure. 
Accordingly, the approximating density functions become increasingly spiky with grow¬ 
ing value of N; see Figure [TT] for an example. Mathematically, we speak of a limit in the 
vague topology. However, the corresponding distribution function F^^\x) = fo f^^Hx)dx 
(which corresponds to the integrated diffraction intensity) converges and possesses a contin¬ 
uous limit; compare the bottom part of Figure [TTJ The limit function can be calculated and 
expressed as an explicit Fourier series; several examples are shown in Figure [T2j 

While this case has no point spectrum (the trivial Bragg peak at 0 being absent due to 
our balanced choice of weights, corresponding to zero average scattering strength), it is by 
no means featureless. In fact, there are peaks that grow with certain scaling exponents (in 
terms of the system size) at certain points in Fourier space (the most prominent examples 
can be found at rational positions | and | in Figure [TT]l . while the growth is not well defined 
at uncountably many other positions (due to the non-convergence of the density functions); 
see Baake, Grimm &: Nilsson (2014) for a detailed account of the classic Thue-Morse case. 

Clearly, the generalised Thue-Morse systems possess hierarchical order, although this is 
not reflected in a pure point component in their diffraction measures. However, this ‘hidden’ 
order is visible in other correlations. Explicitly, it can be revealed by looking at the two- 
point correlations of pairs rather than of single points. Looking at pairs can be described by 
considering the image of the bi-infinite fixed point word v under a sliding block map 4>, which 
maps pairs of adjacent letters to a or 6 according to whether they are equal or different, so 
(j): 11,11 1 -^ a, ll,ll I—)• 6; see Figure [13] for an illustration. 




m 

>|a< 


abmp^ababamaababmbaaaba 


pi 

me 


iiiiiii... 

baaabaa... 


Figure 13: The action of the sliding block map cj) on a Thue-Morse word produces a period 
doubling word. 
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This maps the set of generalised Thue-Morse words to bi-infinite words which are locally 
indistinguishable to fixed point words of the generalised period doubling substitution 

a ^ b^-^ab^-% b ^ b^-^ab^-^a, 

which reduces to the period doubling substitution (with a = 1 and 6 = 0) in the case k = I = 1. 
This map is globally 2:1, meaning that there are precisely two generalised Thue-Morse words 
that are mapped onto the same generalised period doubling word. This is most easily seen by 
noticing that, when going backwards from a generalised period doubling word, there is a single 
free choice for one letter a or 6 at one position, where either preimage can be chosen, after 
which all other preimages are uniquely determined (due to the overlap of adjacent pairs). As 
the generalised period doubling substitution has a coincidence in the sense of Dekking (1978), 
it is pure point diffractive, as discussed above for the (standard) period doubling case. In fact, 
the corresponding point sets are model sets, this time with (k + £)-adic numbers as internal 
space, and the pure point diffraction is supported on the set which contains all inverse 

powers of {k + £) as generating elements. 

In the language of dynamical systems, the dynamical system (where the dynamics is given 
by shifting the sequence by an integer) corresponding to generalised period doubling words 
constitutes a factor of the dynamical system associated to the generalised Thue-Morse words. 
Here, the word factor refers to the fact that it is the image under the sliding block map (j). 
What happens in this case is that the diffraction spectrum of the factor (the generalised period 
doubling system) picks up a non-trivial point spectrum, which is ‘hidden’ in the Thue-Morse 
system, in the sense that it does not show up in its diffraction spectrum (even in the case 
of general weights). However, this pure point spectrum is part of the so-called dynamical 
spectrum of the Thue-Morse system, where the dynamical spectrum refers to the spectrum 
of the operator which generates the translation action; see Queffelec (2010) for details. The 
dynamical spectrum is, in general, richer than the diffraction spectrum. This can be intuitively 
understood because diffraction, as the Fourier transform of the autocorrelation, only measures 
two-point correlations, while the dynamical spectrum ‘knows’ about more general properties 
under the shift action, so effectively can probe higher correlations. We shall come back to 
this point at the end of this section. 

5.3 Order and absolutely continuous diffraction 

Absolutely continuous (‘diffuse’) diffraction is commonly associated with randomness. Indeed, 
stochastic systems typically show absolutely continuous diffraction; the simplest case is the 
Bernoulli shift, based on a random sequence 

A = (...,A_2,A_i,Ao,Ai,A2,...) G {±1}^ 

of independent and identically distributed (i.i.d.) events with outcome ±1, with probability p 
for outcome 1 and probability 1—p for outcome — 1. The Bernoulli shift has (metric) entropy 
H{p) = —plog{p) — (1—p) log (1—p), which is greater than zero except for the deterministic 
limiting cases p = 0 and p = 1. All the examples discussed earlier were deterministic sequences 
with zero entropy. 
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Figure 14: Schematic representation of the Dirac comb of the Rudin-Shapiro chain with 
weights ± 1 , which is again ‘balanced’ in the sense that the average scattering strength is 
zero. 


A random sequence X G {±1}^ leads to a Dirac comb a; = ''^hich is a 

translation bounded random measure with autocorrelation 7 g = Ylmei, % (™) ■ The auto¬ 

correlation coefficients are, almost surely (in the probabilistic sense, so with probability 1 ), 
given by 




1 , m = 0 , 

( 2 p— 1 )^, m 7 ^ 0 . 


as a consequence of the strong law of large numbers. The corresponding diffraction measure 
then, almost surely, is given by 


7^ = (2p - if 6^ + 4p(l -p)X, 

which contains both pure point (for p and absolutely continuous components (for p 0 
{0,1}). The pure point part vanishes when both weights appear with equal probability, while 
the absolutely continuous part vanishes in the two deterministic, periodic cases. 

It is, however, possible to construct deterministic systems with absolutely continuous 
diffraction as well. The paradigm for this situation is the Rudin-Shapiro chain (Rudin 1959, 
Shapiro 1951, Queffelec 2010). Its binary version w G {±1}^ can be defined by initial condi¬ 
tions 'ic(— 1 ) = — 1 , w{0) = 1 , and the recursion 


w{4n + i) 


w{n), for i G {0,1}, 

{—lf~^^w{n), forl'G{2,3}. 


( 21 ) 


A schematic representation of the corresponding Dirac comb is shown in Figure [Ml By 
considering the recursion relation for autocorrelation coefficients induced by Eq. m , in a 
similar way as for the generalised Thue-Morse case above, on can show (Baake & Grimm 2009) 
that the autocorrelation has the simple form 7^3 = <5q, which means that all correlations (apart 
from the trivial case with distance zero) average to zero along the chain. According to the 
two-point correlations, the Rudin-Shapiro chain hence looks completely uncorrelated, exactly 
as the random chain with probability p = ^. As a consequence, the diffraction measure is 
Lebesgue measure, = A, which is clearly absolutely continuous with respect to itself. This 
means that the diffraction intensity is constant in space, and hence completely featureless, 
reflecting the complete absence of two-point correlation in the structure. This example shows 
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that two very different systems such as the p = ^ Bernoulli chain with entropy log(2) and 
the completely deterministic binary Rudin-Shapiro chain (with zero entropy) can produce 
the same autocorrelation and hence the same diffraction measure. Such structures are called 
homometric (Patterson 1944) and show that the inverse problem does not have a unique 
solution in general. 

In fact, the situation is worse than that, as from the results above one can construct an 
entire one-parameter family of stochastic Dirac combs which all are homometric with the 
Rudin-Shapiro chain. This is done by the Bernoullisation procedure (Baake &: Grimm 2009). 
Applying it to the Rudin-Shapiro Dirac comb, the weight at any position is changed randomly 
with probability 1 — p, resulting in 

uj — ^ Wj Xj 5j , 
jez 

with the Rudin-Shapiro sequence w E {±1}^ and the random sequence X E {±1}^ as defined 
above. This is a ‘model of second thoughts’ in the sense that, when starting from a Rudin- 
Shapiro sequence, weights are randomly changed with probability l—p independently at each 
position along the chain. We thus can continuously interpolate between the Rudin-Shapiro 
chain with entropy 0 and the P = \ Bernoulli chain with entropy log(2), with all systems 
sharing the same absolutely continuous diffraction. 

It is interesting to note that the Rudin-Shapiro chain, like the generalised Thue-Morse 
chains above, possesses a ‘hidden’ limit-periodic order that is revealed when looking at an 
appropriate dynamical factor. Using the same sliding block map (p as above, one obtains once 
more a factor with pure point diffraction spectrum, in this case supported on Z[^], as for the 
period doubling case; see Baake &: Grimm (2013) for details. Glearly, this does not happen 
for the stochastic chain, which does not have any ‘hidden’ order. 

5.4 Discrete structures with continuous symmetry 

An interesting (and still somewhat mysterious) class of structures is provided by discrete 
systems which possess a continuous symmetry. The paradigm for such a structure is the 
Gonway-Radin pinwheel tiling (Radin 1994). It is an inflation tiling based on a single trian¬ 
gular prototile of edge lengths 1, 2 and \/5 together with its reflected version. The inflation 
rule is shown in Figure [lU it consists of a linear rescaling by the inflation factor \/5 (first 
step) and the dissection of the inflated triangle into five copies of the original prototile (second 
step), where both orientations occur. The reflected rule applies to the reflected triangle, and 
hence the tiling is reflection symmetric. A realisation of the tiling is shown in Figure [T6l 

What makes this inflation rule special is the rotation it introduces between copies of the 
prototiles. This rotation by an angle i? = — arctan(^) is incommensurate with vr, and as a 
result introduces new, independent directions under inflation. Iteration of the inflation rule 
on an initial patch thus leads to patches comprising an exponentially increasing number of 
triangles occurring in a linearly growing number of independent directions. In the limit of an 
infinite tiling, triangles appear in infinitely many different orientations. While for the familiar 
cases of Penrose-type tilings inflation rules produce tilings with discrete (in the Penrose case 
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decagonal) symmetry (in the sense that the tiling space defined by the fixed point tilings has 
decagonal symmetry; see Baake &: Grimm (2013) for details), the pinwheel inflation produces 
a tiling space with complete circular symmetry (Radin 1994, Radin 1997, Moody, Postnikoff 
(fe Strungaru 2006). As a consequence, its diffraction is circularly symmetric as well, and 
hence cannot have any pure point component apart from the trivial Bragg peak at the origin. 

In fact, the rotation is rather special, because it is a coincidence rotation for the planar 
square lattice, as tan(tl) = — ^ is rational; see Baake (1997) for background. This property 
is behind the observation that the point set of pinwheel reference points can either be seen 
as a subset of rotated square lattices or a subset of scaled square lattices, with scaling by 
inverse powers of 5 (Baake, Frettldh & Grimm 2007a), which makes it possible to draw con¬ 
clusions on the diffraction spectrum by using a radial version of Poisson’s summation formula. 
This provides evidence that the diffraction consists of sharp rings, and that it is surprisingly 
similar to a toy model of powder diffraction of square lattice structures (Baake, Frettloh & 
Grimm 2007a, 2007b). A diffraction measure supported on sharp rings in the plane is sin¬ 
gular continuous, and it is clear that the diffraction of the pinwheel tiling contains singular 
continuous components of this type; however, to date there is no complete characterisation 
of the diffraction spectrum of this example. Results of numerical investigations suggest that 
an absolutely continuous component may also be present. An approximation of the radially 
averaged diffraction is shown in Figure fT71 

While the pinwheel tiling may seem a rather exotic structure, it is generated by a quite 
simple inflation rule with only a single shape up to congruency. There are many other 
structures of this type; see Frettloh (2008) for some examples. 

5.5 Diffraction versus dynamical spectra 

The examples of the Thue-Morse and Rudin-Shapiro systems show that systems can possess 
‘hidden’ order that does not manifest itself by a pure point component in the diffraction 
pattern. However, this order can show up in the dynamical spectrum, which is related to the 
analysis of the translation action on the structure. There is a close relationship between these 
two spectral quantities — indeed, the first proofs of the pure point diffractivity of model sets 
employed the link to dynamical spectra, using the results that the diffraction spectrum is pure 



Figure 15: Inflation rule for the pinwheel tiling. The dot marks the reference point, and the 
shading emphasises that the particular triangle is in the original position and orientation, 
ensuring that repeated application of the inflation rule produces a fixed point tiling. 
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Figure 16: A building at Melbourne’s Federation Square featuring a pinwheel tiling fagade. 

point if and only if the dynamical spectrum is. In general, however, the dynamical spectrum 
can be richer (van Enter & Migkisz 1992), and the Thue-Morse and Rudin-Shapiro systems 
are examples; in both cases, the dynamical spectrum contains the pure point component Z[^] 
which arises because both examples stem from primitive substitutions of constant length 2 
(in the Rudin-Shapiro case, the underlying substitution employs four different letters, and 
the binary system is derived from this by identifying pairs of letters; see Baake &: Grimm 
(2013) for details). 

A particularly simple yet striking example, originally suggested by van Enter, is discussed 
in Baake &: van Enter (2011). It considers the set of certain conhgurations of ±1 on the integer 
lattice Z. The allowed configurations w G {±1}^ are obtained by partitioning the lattice into 
pairs of neighbouring points (there are two ways to do this), and then randomly assigning 
to each pair either the values (+1,-1) or (—1,+1). Turning a configuration w into a signed 
Dirac comb with weights Wi G {=tl}, it is easy to show, by an application of the strong law 
of large numbers, that the autocorrelation is (almost surely) given by 7 = Jq — ^(5;^ + (5_i). 
The corresponding diffraction measure is then 

7 = (1 — cos(27r/c)) A , 

and hence purely absolutely continuous, where the Radon-Nikodym density relative to A 
is written as a function of k. However, the dynamical spectrum of this system contains 
eigenvalues (hence a pure point part), reflecting the order in the system imposed by the 
‘dimer’ condition on pairs. This can be revealed by considering a block map similar to to the 
map (p used above. Explicitly, setting Ui = —tCjtCj+i for f G Z maps tc to a new sequence u, 
which (almost surely) has the diffraction measure 
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Figure 17; Approximation of radial diffraction intensity I{k) for the pinwheel diffraction 
(black line), based on data from a finite system. The grey columns indicate the sharp rings 
observed in a toy model of powder diffraction from a planar square-lattice structure, with the 
relative scale adjusted according to the first main peak; see Baake, Frettloh & Grimm (2007a, 
2007b) for details. 

see Baake &: van Enter (2011) for details. The ‘dimer’ structure is reflected in the presence 
of the pure point part supported on ^Z, which also is the entire point part of the dynamical 
spectrum. 

This example again shows that the ‘hidden’ order can also be seen in diffraction, but not in 
the original system. Note that simply changing the weights of the scatterers will not achieve 
this, although it may contribute a trivial Bragg part. However, choosing a suitable factor (or 
a family of factors) as an image of a continuous map such as the sliding block map (p used 
above, makes it possible to detect the ‘hidden’ order via its diffraction. That this is a general 
property of the relation between dynamical and diffraction spectrum is a recent non-trivial 
insight; see Baake, Lenz & van Enter (2013) for the latest developments in this direction. 

6 Conclusions 

The discoveries of incommensurately modulated and aperiodically ordered solids in the twenti¬ 
eth century (de Wolff 1974, Janner &: Janssen 1977, Shechtman, Blech, Gratias & Gahn 1984, 
Ishimasa, Nissen &: Eukano 1985) have changed our view of crystallography. Grystallography 
is no longer restricted to the analysis of lattice periodic arrangements of atoms or molecules, 
but takes a broader view which includes certain aperiodically ordered structures, such as 
incommensurate crystals and quasicrystals. The definition of a crystal has been amended to 
reflect this broader view. 

The definition of a crystal is based on the currently known catalogue of periodic and 
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aperiodic crystals. We presently do not know of any materials that have aperiodically or¬ 
dered structures beyond incommensurate crystals (including composite structures) and qua¬ 
sicrystals. For the latter, so far only symmetries corresponding to the smallest embedding 
dimension (in the sense of model sets) have been observed, with octogonal, decagonal and 
dodecagonal quasicrystal planes corresponding to projections from four-dimensional periodic 
structures, and icosahedral quasicrystals being described by projection from six-dimensional 
lattices. However, there is no a priori reason that excludes other symmetries completely, or 
indeed aperiodically ordered structures that are not described by model sets obtained from 
projections of a lattice in a finite-dimensional Euclidean space. 

The definition of a crystal also reflects the current lack of understanding of what consti¬ 
tutes order in matter (and more generally), and in this sense should be seen as a working 
definition that may well need to be revised in the future. In crystallography, order is linked to 
diffraction, which makes sense because diffraction is the method of choice to experimentally 
determine the structure of a solid. The examples discussed above demonstrate that there are 
ordered structures which are not captured by the current definition, either because their pure 
point diffraction fails to be finitely generated, or because they do not have any non-trivial 
point component in their diffraction. While we do not know whether such structures are 
realised in nature, it should become possible to manufacture such materials with purpose- 
designed structures and properties. In this sense, these structures are relevant and should be 
considered to be within the realm of crystallography. 

From a mathematical point of view, a more satisfying attempt at defining order might 
employ the dynamical spectrum, which is a generalisation of the diffraction spectrum. The 
results above are in line with the intuition by van Enter & Mi§kisz (1992) that an appar¬ 
ent disorder at an ‘atomic’ scale could be accompanied by order at a ‘molecular’ scale, with 
diffraction of derived factor structures probing the latter. While diffraction itself only mea¬ 
sures the averaged two-point correlations in a structure, the dynamical spectrum probes the 
repetitivity of a structure under translations, and hence also higher-order correlations, which 
generally can distinguish homometric systems (Griinbaum & Moore 1995). While these are 
not necessarily directly accessible by experiment, the additional information contained in the 
dynamical spectrum is, in principle, encoded in diffraction spectra of derived systems; see 
Baake, Lenz & van Enter (2013) for recent developments on establishing this connection. 
Defining order via a non-trivial pure point component of the dynamical spectrum would in¬ 
clude structures such as the Thue-Morse and Rudin-Shapiro systems, though presumably 
examples of pinwheel type (for which the dynamical spectrum is not known) would be ex¬ 
cluded. In this sense, it probably is still not completely satisfactory to capture all possible 
manifestations of order, but it may provide a first step towards a better understanding. 

In this paper, the discussion was limited to deterministic systems, apart from the brief 
excursion on the Bernoulli chain. Clearly, moving to partially ordered systems, which contain 
an element of stochastic disorder, is relevant as well. Not only does even the most perfect 
crystal contain some amount of disorder, but there are also entropically stabilised structures 
with intrinsic configurational disorder, among them many quasicrystalline phases. In this 
context, the notion of ‘entropic order’ is relevant, which has been investigated in statistical 
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physics, in particular with respect to the physics of glasses; see, for instance, Kurchan & 
Levine (2011), Sasa (2012a, 2012b) and Wolff &: Levine (2014) for recent work along these 
lines. 

Although the importance of random tiling structures was pointed out early on (Elser 
1985), and while there is some good heuristic information from scaling considerations (Henley 
1999), there are as yet very few mathematically rigorous results for non-trivial random tiling 
structures in two or more dimensions, the only examples known being related to solvable 
models of lattice statistical mechanics (Baake &: Hoffe 2000). We refer to Baake, Birkner & 
Grimm (2015) for a recent review on what is known about the diffraction of partially ordered 
and stochastic systems. 
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